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We study correlation effects on the transport through a quantum dot superlattice 
using a two-dimensional Hubbard model connected to two noninteracting leads. To 
calculate the zero-temperature conductance away from half-filling, we have used the 
non-magnetic solution of the Hartree-Fock approximation for the unperturbed Green's 
function, and then take into account the electron correlation beyond the mean field 
theory through the second-order self-energy corrections with respect to the residual 
interaction. When the value of the onsite potential ea or the repulsion U is changed, 
the conductance shows a number of peaks corresponding to the resonance states passing 
through the Fermi energy. The behavior of the resonance states for correlated electrons 
can be investigated using the eigenvalue of the effective Hamiltonian for a quasi-particle 
of a local Fermi liquid. We calculate the eigenvalue using the second-order self-energy, 
and show explicitly that it can be compared with the conductance peak. 
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1. Introduction 

Systems of nanometer scale and artificial structures, such as quantum dots and quantum wires, 
have been a subject of current interest. In these systems the Coulomb interaction seems to play an 
important role on the transport properties. For instance, the Kondo effect in quantum dots has been 
studied intensively 1-5 ) and the experiments have been shown to be in qualitative agreement with 
theoretical predictions. 6, 7 ) The artificial molecule which can be realized by arranging two or more 
quantum dots is another example. 8, 9 ) Also, an interesting quantum dot superlattice showing the 
ferromagnetism or c?-wave superconductivity has been designed theoretically. 10, n ) The proposed 
superlattice consists of a two-dimensional network of quantum wires, and is described by a Hubbard 
model that is derived based on the local density approximation. 

In theoretical studies of the transport properties of strongly correlated electrons, accurate nu- 
merical approaches such as the numerical renormalization group 12,13 - 1 and quantum Monte Carlo 
methods 14, 15 ^ have been used successfully. However, the analytic methods, which are complemen- 
tary to the numerical ones, are also required to clarify the underlying physics comprehensively. 
The perturbation theory in the Coulomb interaction is one of such approaches. For instance in 
the Kondo problem, 16 ' the perturbation expansion with respect to U has shown to work quite 
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well to describe the low-energy Fermi-liquid behavior. 17,18 ) The perturbation theory can also de- 
scribe properly the low-voltage behavior of the out-of-equilibrium Kondo system under finite bias 
voltages. 5 ' 19 ^ Furthermore, based on the analytic properties of the vertex function, the transport 
equation for the Fermi liquid derived by Eliashberg 20 ) has been shown to have a link to the Landauer 
formula for the conductance of the correlated electrons. 21 ) 

We have been studying the transport properties of a quantum dot superlattice based on a Hub- 
bard model coupled to two noninteracting leads. 22 ) A schematic picture of the model is illustrated 
in Fig. 1: the system consists of a Hubbard model of Ac (= A x M) sites and two semi-infinite 
noninteracting leads of M channels at left and right, where N and M are the size in the x- and 
y-direction, respectively. We have considered the square lattice for simplicity. At this moment this 
model itself might have no direct correspondence to experiments. Nevertheless, it seems to contain 
essential features of two-dimensional network of quantum dots or atomic systems, and we believe 
that it could be realized some day in near future. In this kind of the system that coupled to the 
noninteracting leads, the contribution of the vertex corrections on the conductance vanishes at zero 
temperature. 23 ) Thus, at T = 0, the conductance can be written in terms of the single-particle 
Green's function at the Fermi energy oj = 0. In previous work we have calculated the conductance 
at half-filling, where the average number of the electron per site is one, using the order U 2 self- 
energy. The conductance shows some peaks as a function of U when the hopping matrix elements 
are anisotropic. We found at half-filling that the position of the resonance peaks can be analyzed 
through the eigenvalues of the effective Hamiltonian for free quasi-particles. However, in actual 
quantum dots the half-filled situation is expected only for a special value of the gate voltage, which 
can be related to the onsite energy of the theoretical model. Therefore, real systems are usually 
in away from half-filling, so that it should be examined on the same basis. 

The purpose of this paper is to study the correlation effects on the transport properties away 
from half-filling. In this case, the average number of electrons per site is not unity, and the spatial 
distribution of the charge is not uniform. Therefore, as the density functional theory states, 24 ) 
the charge distribution should be determined accurately to obtain the ground-state properties cor- 
rectly. Note that the charge distribution is uniform at half-filling, so that we did not have to take 
any care over the details of the distribution in the previous study. In the present study, we have 
used a non-magnetic solution of the Hartree-Fock approximation to determine the unperturbed 
Green's function, and then take into account the correlation effects beyond the mean field approx- 
imation through the second-order self-energy corrections with respect to the residual interaction. 
Near half-filling, the result of the conductance for even Ac shows a valley-like structure in the 
ed dependence. Comparing with the Hartree-Fock results, we see that the self-energy corrections 
make the valley deep and wide. This is consistent with the tendency towards the development of 
a Mott-Hubbard insulator gap, and we have examined this feature both for one-dimensional chain 
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and two dimensional lattice. For odd Nq the correlation effects make the conductance peak, which 
is seen in the dependence near half- filling, flat. This can be regarded as a manifestation of the 
Kondo effect, and we have carried out the calculations for odd Nq using a one-dimensional chain. 

In the finite system connected to reservoirs, the conductance shows a peak when a resonance 
state pass through the Fermi energy. Especially, away from half-filling the peak structure can be 
seen by changing the value of e^. The behavior of the resonance states can also be analyzed in 
terms of an extended version of a local Fermi-liquid theory. The presence of the noninteracting 
leads makes it possible to construct an effective Hamiltonian of a bilinear form, which reproduces 
exactly the conductance and the total charge displacement of the interacting system at T = 0. In 
this paper we demonstrate using the order U 2 self-energy corrections that the eigenvalue of the 
effective Hamiltonian can be used also away from half-filling for investigating the position of the 
resonance states. 

Horvatic and Zlatic have shown for the single impurity Anderson model that the radius of the 
convergence of the perturbation expansion with respect to U decreases rapidly when the system 
goes away from half-filling. 25 ) The same thing happens also for systems consisting of a number of 
interacting sites, and away from half-filling the second-order perturbation theory is applicable for 
relatively small values of U compared to that at half-filling. Therefore, in the present work we have 
concentrated on the weak interaction cases, and have checked the applicability by comparing the 
occupation number which are obtained with the two different ways; i) integration of the density 
of states up to the Fermi level, and ii) the Friedel sum rule. The results obtained with these two 
methods agree well for small U. 

In §2, we describe the model and formulations for studying the correlation effects away from 
half-filling. In §3, we present the results of the conductance and the average number of electrons, 
and discuss the interpretation based on the effective Hamiltonian. Summary is given in §4. 

2. Model and Formulation for away from Half-filling 

In the system without the translational invariance, the spatial distribution of the charge is 
usually not uniform, except for the electron-hole symmetric case where the average number of elec- 
trons on each site is fixed as one. Therefore, away from half-filling, the charge distribution depends 
on the parameters of the Hamiltonian, and an extra care is necessary to construct the unperturbed 
part appropriately. In this section, we describe the outline of the approach: the unperturbed part 
is determined with the Hartree-Fock approximation assuming that the ground state to be param- 
agnetic, and then the electron correlation is taken into account using the perturbation expansion 
with respect to the residual interaction. 

Fig. 1 shows the schematic picture of the model, where the lattice sites (quantum dots) are 
labeled with the coordinate r = The interacting region, where the onsite repulsion U is finite, 
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consists of M rows in the y-direction and N columns in the x-direction: t x and t y are the nearest- 
neighbor hopping matrix elements in the x- and y-directions, respectively. At the columns of x = 1 
and x = N, the interacting region is connected to noninteracting leads with the mixing matrix 
elements vi, and vr, which are assumed to be uniform in the y-direction. The total Hamiltonian H 
is given by 
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Here c| creates an electron with spin a at r = and nj Ji0 - = c| ■ CT Cj • CT . The periodic 

boundary condition is assumed along the y-direction c^m+i.o- = Cj,i,<r- In eq. (5), is the bare onsite 
energy, and (wjj)o = So-( n i,:/>)o is the non-magnetic solution of the Hartree-Fock approximation: 
i n i,j,])o = ( n i,j, i)o, this average generally depends on the position The Hamiltonian for 

the two leads, Hl and Hr, are written in terms of the one-particle energy e Xk and corresponding 
eigenfunction (pxk( r x) f° r A = L, R. The operator a X j a = J2k a \ka ( i ) Xk( r x j) appearing in the 
mixing Hamiltonian 7Y m i x annihilates an electron at the interface r x -. We will be using units h = 1 
unless otherwise noted. 

The single-particle Green's function is defined by 
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Here [3 = 1/T, e n = (2n + l)ir//3, c rcr (r) = e rW c ro .e rW with r = and (• • • ) denotes the 

thermal average Tr [ e~^ w • • • ] /Tr e _/3W . Since J7 is acting only in the central region (C), the Dyson 
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G rr i(ie n ) = G° rr ,(\e n ) + ^ G° rri (i£ n ) £ rir2 (ie n ) G r2r /(ie n ) . 

ri,r 2 ec 



(9) 



Here G^ r ,(ie n ) is the unperturbed Green's function corresponding to 7i°, which is determined with 
the non-magnetic solution of the Hartree-Fock approximation, so that 

2 r° 
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The value of (rar)o is determined self-consistently from eqs. (5) and (10). To specify the various 
types of the Green's functions, we use the symbol + (— ) in the superscript as a label for the 
retarded (advanced) function, i.e., = G rr i(uj ± i0 + ). In eq. (9), ^ rj ,t is the self-energy 

correction beyond the mean-field theory, and the summations with respect to r\ and r 2 run over 
the Nq = N x M sites in the central region. Therefore, the Dyson equation can be written in a 
matrix form, 
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In eq. (13), 1 is the M x M unit matrix, and h° is a tridiagonal M x M matrix consisting of Eij 
at the j-th main diagonal part and — t y at sub- and superdiagonal parts. In eq. (15), £jj'(z) is a 
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M x M matrix whose (j,f) element is the self-energy £ ry ./ for r = r' = (i',j f )- In eq. (14), 

F\{z) = {F\ : jj>(z)} is the Green's function at the interface of the isolated lead at A = L, R; 
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From this and the hybridization v x , the level width due to the mixing with the leads is given 
by T\(uj) = —v x lmF x (uj) for A = L and R. In the present study, we assume a diagonal form 
F^(oj) = F^(uj) = —in pi taking the local density of states at the interfaces p to be a constant 
independent of u, so that the effects of the mixing can be parameterized by T\ = irpv\. Note that, 
because of the periodic boundary condition in the y-direction, the M x M matrix functions can be 
diagonalized as = Ylm=i Xm^tf'xL Here x m i s the eigenvector for the subbands, and the 
corresponding eigenvalue is given by e m = —2t y cos(27rm/M) for m = 1, 2, . . . , M. 

At zero temperature, the dc conductance can be calculated from the Green's function as 
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where Gu>(z) is a M x M matrix that corresponds to the (i,i r ) partitioned part of Q{z). Fur- 
thermore, at T = 0, the total charge displacement AiVtot can be calculated from the Friedel sum 
rule, 22 ' 28) as 
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Specifically, in the case of the constant density of states p, the Anderson compensation theorem 
holds, 29 ) and then the charge displacement of whole the system coincides with the number of 
electrons in the central region: 22 ) 

N M 

AiVtot = £E< n «>- (20) 

i=l j=l 

These expressions, eqs. (17) and (18), show that at T = the conductance and charge displace- 
ment can be determined by the value of the Green's function at the Fermi energy u = 0. Therefore, 
one can introduce a noninteracting Hamiltonian that leads exactly the same g and AiV to t of the 
interacting system. 23 ) This is due to the property ImS + (0) = at T = 0, and in the central region 
such an effective Hamiltonian can be constructed as 
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H c + ReS (0) . 



(21) 
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Note that the self-energy is defined with respect to the whole system, so that ReS + (0) depends 

not only on U but also on T L and T R through the unperturbed Green's function Q (z). We see in 

^ eff 

the next section that the eigenvalues of 7i c correspond approximately to the peaks of the resonant 
states. 

3. the Second Order Perturbation Theory 

In this section, we calculate the self-energy ReX! + (0) up to the second order with respect to 
Ti 1 ^ defined by eq. (6), and show the results for the conductance and charge displacement. In 
order to determine the unperturbed Green's function G , with the Hartree-Fock approximation 
assuming the paramagnetic ground state, we calculate first the value of Eij by solving eqs. (5) and 

(10) for given and U. Then, we take into account the lowest order self-energy correction due to 
the fluctuation Hq*. It is described by the diagram shown in Fig. 2, and its value at T = and 
to = is given by 

Kr>^ =- u2 rjz w G ° rr,{[£) G ° rr,{[£,) + i£ ° • (22) 

The explicit form of G^, can be derived from eq. (12). Note that the retarded function at u = 
and T = can be obtained from the Matsubara function, i.e., ,(0) = S rT ,/(ie)| e _ > o+- From 
eq. (22), one can also confirm explicitly the property Im£+ ,(0) = 0, at oj = and T = 0. We 
calculate all the Nq x Ac elements of ReS^ 7 (0) carrying out the integrations numerically. Then, 
we obtain the full Green's function G + ,(0) substituting Re£ + ,(0) into the Dyson equation eq. 

(11) . 

In the previous work we have used this method to study the ground state properties at half- 
filling, 22 ) where the spatial distribution of the charge is uniform independent of the details of the 
parameters. In this case the perturbation theory works rather well. 23 ) Specifically, for the single 
Anderson impurity model, which corresponds to a special limit A^c = 1, the order U 2 self-energy 
describes the low-energy Fermi-liquid behavior properly 17 ) and gives the spectral function that 
is qualitatively correct even for rather large U. 16 ^ Away from half-filling, however, the spatial 
distribution of the charge is inhomogeneous, and it makes the applicability of the second-order 
perturbation theory rather worse. As the density functional theory 24 ) teaches us, in order to 
get the ground state correctly, one needs the method that can reproduce the charge distribution 
accurately. We know the correct charge distribution at half-filling but we do not away from half- 
filling. Therefore, away from half-filling the charge distribution is the quantity that should be 
calculated. This makes the situation away from half-filling different from that at half-filling. For 
the Anderson impurity, the second-order perturbation theory has been applied to away from half- 
filling, 25 ) and a modified method has also been examined by several authors. 26,27 ) Horvatic and 
Zlatic has shown for the Anderson model that the radius of the convergence of the perturbation 
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expansion with respect to U decreases rapidly when the system goes away from half-filling. 25 ) We 



see that there is a similar tendency also for Ac > 1 in the following subsection. To check the 
applicability of the approximation, we calculate the total charge in the interacting region with two 
different methods. The first one is using the integration of the local density of states (DOS) up to 
the Fermi level, 



where Q (uj) should be calculated using the order U 2 self energy. The second method is using 
the Friedel sum rule eqs. (18) and (19), and thus the charge is determined by the value at the 
Fermi energy u = 0. The two results agree well for small U, and it gives some insights into the 
applicability. 

In the following subsections, we apply the formulation to the Anderson impurity in §3.1 in order 
to demonstrate how it works in the simplest case. Then, we examine the one-dimensional chain 
in §3.2 and the two-dimensional lattice in §3.3. Specifically, we concentrate on the case Tl = Tr 
(= r), where the system has the inversion symmetry. 

3.1 Anderson impurity 

We consider here the single Anderson impurity that corresponds to the case of Nq = 1 (M = 1, 
N = 1). In this subsection, we use the width of the resonance peak for the noninteracting case as 
the unit of the energy, i.e., A = 2r. 

In Fig. 3, the average number of electrons at the impurity site (n^) is shown as a function of 
U for e^/A = —0.3 and —1.0. The system has the electron-hole symmetry at the value of U where 
(rid) is equal to 1.0. For U / (ttA) < 1.0, the results obtained from the Friedel sum rule (•) agree well 
with that of the integration of DOS (A). The range of U where we have good agreement depends 
on ed, but the dependence is weak for small e^. When the onsite energy is unchanged, the local 
charge (n^) must be a decreasing function of U and converges to a finite constant in the U — > oo. 
Therefore, the upturn seen in Fig. 3 indicates the breakdown of the approximation for large U. In 
Fig. 3, we have also shown the Hartree-Fock results (dashed lines) for e^/A = —0.3 and —1.0. For 
U < — led (or equivalently (rid) > 1-0), the Hartree-Fock result is larger slightly than that of the 
second-order perturbation result, although it might be hard to see this feature in the figure. In 
contrast, for U > — 2e^ (or equivalently (rid) < 1-0), the Hartree-Fock results is smaller than that 
of the perturbation results, and the difference increases with U. These results are consistent with 
that in the seminal work by Horvatic and Zlatic, 25 ) and suggest that the method is applicable for 



In Fig. 4, the result of the conductance is shown as a function of [// (7rA) for several values of 
e^/A = —0.3, —0.7, —1.0, and —1.5. The conductance shows a maximum at half-filling (n^) = 1, 




(23) 



EV(ttA) < 1.0. 
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which corresponds to U = — 2ed- This is because at half-filling the resonance level is situated 
exactly at the Fermi level. The change of the position of the resonance state can be traced through 

•~ eff 

the eigenvalue of the effective Hamiltonian H c . For the Anderson model the eigenvalue is given 
by Ed = ed + U (nd)o/2 + ReS + (0), which is plotted as a function of U in Fig. 4 (c). The perfect 
transmission occurs, when the eigenvalue crosses the Fermi level at u = 0. In the case of the 
Anderson impurity, corresponds exactly to the position of the resonance state. However, for 
larger systems consisting of a number of the interacting sites Ac > 1> the precise position of the 

^ eff 

resonant peaks deviates slightly from the eigenvalue of 7i c because of the additional effect of the 
hybridization V m ix appearing in the right-hand side of eq. (12). 

3.2 One-dimensional chain 

We next consider one-dimensional chains connected to noninteracting leads, which is the case of 
M = 1 and N > 2. The system can be regarded as a model for a series of quantum dots or the atomic 
chain of nanometer size. In the electron-hole symmetric case the second-order perturbation theory 
described above has already been applied to this model. 23 ) The results obtained at half-filling show 
the even-odd oscillatory behavior: the perfect transmission due to resonant tunneling occurs for 
odd N, while for even N the conductance decreases with increasing N showing a tendency towards 
the development of a Mott-Hubbard insulator gap. In this subsection, we study the ground state 
properties away from half-filling. In the following, we use the hopping matrix element t, instead 
of A, as the energy scale, and take the mixing parameter to be T/t = 0.75 in the numerical 
calculations. 

To check the applicability of the perturbation theory in the one-dimensional case, we have 
compared the number of electrons in the interacting region AA tot with the Friedel sum rule and 
the integration of DOS up to the Fermi level as we have done for the Anderson impurity in the 
above. Specifically, we have examined the chain of length N = 3 and N = 4 taking the onsite 
energy to be ea/t = —0.3 and —1.0. Although the results are not shown here, we have obtained 
the figures, the feature of which is quite similar to that of Fig. 3. The results of the two different 
methods agree well at U/(2irt) < 0.7 both for N = 3 and N = 4. Therefore, in this subsection we 
discuss the results obtained in the small U region; U/(2irt) < 0.7. 

In Figs. 5 and 6, the ground state properties of the chains of length N = 3 and N = 4 are 
shown, respectively, where the onsite energy e^/t is taken to be (•) —0.3, (o) —0.7, and (a) —1.0. 
The system of N = 3 and 4 can be regarded as the simplest examples of the odd and even chains, 
respectively, which capture the essence of the one-dimensional network. Note that the system has 
the electron-hole symmetry when U + 2 = 0. It corresponds to the value of U, at which the local 
charge becomes AJV tot = N and the conductance shows a maximum (minimum) for odd (even) N. 
When an eigenvalue in the middle, which is shown in the panel (c) of Figs. 5 and 6, crosses the Fermi 
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level, the conductance shows a peak. For N = 3, this happens only at half-filling in the situations 
shown in Fig. 5 (c), and the position of the conductance peak coincides exactly with the intersection 
of the central eigenvalue and the Fermi level. In the case of N = 4, the conductance shows some 
peaks away from half-filling. We see two peaks in the curves for (o) e^/i = —0.7 and (A) —1.0 in 
Fig. 6 (a). These peaks in each of the curves correspond to the eigenvalues which cross the Fermi 
level at u = successively in Fig. 6 (c). As mentioned in the above, the position of the resonance 
peaks deviates from that of the eigenvalue because of the energy-shift due to the hybridization. 
The conductance for (•) e^/t = —0.3 in Fig. 6 (c) shows a single peak because only one eigenvalue 
among the four crosses the Fermi level. In the noninteracting case this feature can be seen also in 
the ui dependence of the transmission probability, which is shown in Fig. 7 for e^/i = —0.3. In the 
solid line (for N = 4) two peaks among the four are situated below the Fermi level u = 0, and they 
move towards the Fermi level when U is switched on. 21 ) These examples show that the eigenvalues 
contain information about resonant states situated near the Fermi level. We note that away from 
half-filling not only the off-diagonal part 23 ) but the diagonal part of ,(0) is finite and plays an 
important role on the charge distribution. Furthermore, we note that an analogous description has 
been used in the translational invariant systems, for which the self-energy can be diagonalized in 
the k (wave- vector) space and the dispersion for the quasi-particles has a correction through the k 
dependence of the self-energy. 

So far, we have discussed mainly the U dependence of the conductance and other quantities. 
In real quantum dot systems, however, the onsite energy can be controlled by the gate voltage. 
Therefore, the dependence should also be examined. In Fig. 8 the conductance for (a) N = 3 
and (b) iV = 4 is shown as a function of taking the Coulomb interaction to be U/(2irt) = 0.5. 
The conductance is symmetric with respect to the half-filled point = — U/2 (~ — 1.6 1) because 
of the electron-hole symmetry of the model. To see the effects of the self-energy correction due to 
the residual interaction, in the figures we have also shown the Hartree-Fock results (o) obtained 
using the unperturbed Green's function. In the case of N = 3, the conductance shows a peak at 
half-filling, and the self-energy correction makes it broader. This is due to the Kondo effect, and is 
consistent with the numerical renormalization group results 13 -* for the Anderson impurity. On the 
other hand, in the case of A = 4, the self-energy correction makes the valley of the conductance 
around the half-filling deep and wide. This can be explained as a tendency towards the development 
of a Mott-Hubbard insulator gap. 

3.3 Two-dimensional lattice 

In this subsection we consider the two-dimensional Hubbard model, as a model for a quantum 
dot superlattice. In the numerical calculations, we take the hopping matrix elements to be isotropic 
t x = t y (= t), and take the mixing parameter to be T/t = 0.75. Specifically, we examine a system 
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of the size M = 4, N = 8: the number of the interacting sites is Nq = 4x8. 

In Fig. 9, the conductance and the number of electrons in the interacting region are shown for 
several values of ed/t = —0.3, —0.7, and —1.0. In the lower panel (b) we see that the local charge 
AJVtot, which is obtained using the Friedel sum rule, behaves properly as a decreasing function 
of U at J7/(27r£) < 0.9. In this range of U, the results obtained with the integration of the DOS 
(not shown) agrees with the ones in the panel (b), and we can conclude that the second-order 
perturbation theory is applicable. The conductance shows the peak when a resonance peak passes 
through the Fermi level, and shows a minimum at half-filling, i.e., at U = —2ea or AA^ot = 32. Due 
to the periodic boundary condition in the y-direction, we can see the contribution of each subband 
separately. In the case of M = 4, the subband energy e m = —2t y cos(27rra/M) is given by ei = 2t y , 
en = em = 0, and erv = — 2t y , where the labels are assigned as I (m = 2), II (m = 1), III (m = —1), 
and IV (m = 0). In Fig. 10 the results of the conductance are shown separately, where ea/t is taken 
to be (a) —0.3, (b) —0.7, and (c) —1.0. The resonance states belonging to the mode I (•) and 
mode IV (A) are sharp, while the peaks in the curve for the mode II and III (o) are broad. Note 
that the two degenerate subbands, II and III, make the same contribution. Comparing the panels 
(a)-(c), we see that the number of peaks increases when the onsite energy moves deep inside 
the Fermi level. This is because the resonant states below the Fermi level can cross the Fermi level 
when the repulsive interaction U increases. In Fig. 11, the local charge is plotted separately for 
each subband, where ea/t is taken to be (•) —0.3, (o) —0.7, and (A) —1.0. Naturally, the number 
of electrons in the subbands at low (high) energies is large (small). A step-like behavior is seen for 
the mode I and IV. This is because the occupation of the sharp resonance state changes rapidly 
from 2 to 0, when it crosses the Fermi level. 

The feature of the resonance peaks can be compared with the transmission probability in the 
noninteracting case U = which is shown in Fig. 12 for ea/t = —0.3. In this figure, only the 
contribution of the mode I (solid line) and mode IV (dashed line) is plotted. The contribution of 
the mode II and III, which is not shown in Fig. 12, is identical to the dashed line in Fig. 7 because 
of en = em = 0. The sharp peaks seen in the curves for the mode I and IV in Fig. 10 correspond to 
the ones at the band edge of the noninteracting subband-spectrum. Furthermore, the broad peaks 
seen in Fig. 10 for the mode II and III correspond to the ones in the middle of the noninteracting 
subband-spectrum. Therefore, when the interaction is switched on, the resonance peaks which 
already exist in the noninteracting system still remain changing their position and shape gradually. 

^ eff 

As discussed in the previous subsections, the eigenvalue of the effective Hamiltonian 7i c defined 
by eq. (21) has an additional information about the position of the resonance peaks. In Fig. 13, 
the U dependence of the eigenvalues is shown for e^/t = —1.0. Among N (= 8) eigenvalues in each 
of the modes, those staying near the Fermi level contribute to the conductance. Comparing with 
the conductance shown in Fig. 10 (c), the resonant peaks seen in the contribution of the mode I 
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(IV) corresponds to the lowest (highest) two eigenvalues which cross the Fermi level in Fig. 13 (a). 
Also, the broad peaks seen in the contribution of the mode II and III in Fig. 10 (c) correspond to 
the three intersections at Fermi level in Fig. 13 (b). Therefore, from the eigenvalue just below the 
Fermi level, one can anticipate when the next peak will contribute to the tunneling current. This 
would not be possible if one sees only the conductance. 

In Fig. 14 the conductance of the same two-dimensional lattice is shown as a function of the 
onsite energy taking the Coulomb interaction to be U/(2irt) = 0.5. Note that, as mentioned 
in the previous subsection, can be controlled by the gate voltage. As seen in the lower panel 
(b), the contributions of the mode I and IV are symmetric with respect to the point = — U/2 
(~ — 1.6 1) which corresponds to the half-filling. Around this point an energy gap of the order of U 
opens up between these two subbands, which is not present basically in the noninteracting case in 
Fig. 12. However, since the mode II and III carry the current even in this energy gap region, the 
total conductance is finite in the plotted range of in the upper panel (a). For comparison, the 
Hartree-Fock result of the total conductance is also shown (dashed line) in Fig. 14 (a). Due to the 
self-energy correction beyond the Hartree-Fock approximation, the valley of the total conductance 
around the half-filling becomes deep and wide. Comparing with the one-dimensional chain of 
N = 1 x 4 [see Fig. 8 (b)], this tendency towards the development of a Mott-Hubbard insulator gap 
becomes clear in the two-dimensional lattice of Ac = 4x8. 

4. Summary 

In summary, we have studied the effects of electron correlation on the transport through a 
two-dimensional Hubbard model of finite size connected to two noninteracting leads, as a model for 
a quantum dot superlattice or materials on a nanometer scale. Specifically, we have examined the 
ground-state properties away form half-filling by extending the previous work done at half-filling. 
In the calculations, we determine first the unperturbed Green's function using the non-magnetic 
solution of the Hartree-Fock approximation, and then take into account the second-order self- 
energy corrections with respect to the residual interaction that is not included in the mean field 
approximation. This method works well for small U. However, as it is known for the single 
Anderson impurity, the radius of the convergence of the expansion with respect to U decreases 
rapidly when the system goes away half-filling. 25 ) We have checked carefully the applicability of 
the method through the U dependence of the local charge in the interacting region AAtot, and 
study the transport property at small U where AA tot behaves properly as a decreasing function of 
U. 

We have calculated the zero-temperature conductance through the one- and two-dimensional 
Hubbard systems consisting of several tens of interacting sites. The conductance shows peaks 
caused by the resonant tunneling as a function of the onsite energy which can be controlled by 
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the gate voltage, and as a function of the Coulomb interaction U. Our numerical results show that 
the resonant peaks seen in the conductance have direct correspondence to the eigenvalues of the 
effective Hamiltonian which describes the quasi-particles of the local Fermi liquid. Furthermore, 
from the eigenvalue we can get information about the resonance states inside the Fermi level. 

The formulation used in the present study can be applied to various systems. We have also 
been studying how the randomness in the sample affects the transport through small interacting 
systems: real artificially systems seem more or less to be disordered. We have reported preliminary 
results obtained at half-filling, 30 ) and we are going to study the effects of the randomness away 
from half-filling. 
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M = 1,N = 1 

i i i i 

1 2 

U / (7iA) 

Fig. 3. The number of electrons at the impurity site {rid) of the Anderson model (M = 1, N = 1), where e<j/A is 
taken to be —0.3 and —1.0. The results are obtained with the two different methods which use (•) the Friedel sum 
rule and (A) the integration of DOS up to the Fermi level. The two dashed lines are the results obtained with the 
Hartree-Fock approximation for the same parameters: e^/A — —0.3 and —1.0. 



Effects of Electron Correlation on the Transport through a Quantum Dot Superlattice 



17 




Fig. 4. Ground state properties of the single Anderson model: (a) conductance, (b) local charge (rid), and (c) 
^ e ff 

E d = td + U (n d )o/2 + Re£ + (0) which corresponds to the eigenvalue of "H c . Here the onsite energy e d /A is taken 
to be (•) -0.3, (o) -0.7, (A) -1.0, and (A) -1.5. 
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U / (2nt) 



Fig. 5. Ground state properties of a one-dimensional Hubbard model of odd number sites (M = 1, N = 3): (a) 

conductance, (b) local charge in the interacting region Y2iLi ( n i,j)> an< i ( c ) eigenvalues of W c . Here e<j/t is 

taken to be (•) -0.3, (o) -0.7, and (a) -1.0. 
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U / (27Ct) 



Fig. 6. Ground state properties of a one-dimensional Hubbard model of even number sites (JVf = 1, N = 4): (a) 

conductance, (b) local charge in the interacting region Y2iLi S^li { n i,j)> an d ( c ) eigenvalues of He . Here e<j/t is 
taken to be (•) -0.3, (o) -0.7, and (a) -1.0. 
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co/t 



Fig. 7. Transmission probability vs u for noninteracting chains (M — 1), where ea/t = —0.3 and T/t — 0.75. The 
solid line is for the system of length N = 4, and the dashed line is for N = 8. 
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Fig. 8. The td dependence of the conductance through the one-dimensional chains (M = 1) of length (a) N = 3, and 
(b) N — 4. The solid circles (•) denote the results obtained using the second-order self-energy corrections, and the 
open circles (o) denotes the Hartree-Fock results. The Coulomb interaction is taken to be U/(2irt) — 0.5, so that 
the electron-hole symmetry holds at ej, = —U/2 (~ — 1.6 1). 
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Fig. 9. Ground state properties of the two-dimensional Hubbard model of the size M = 4, N = 8, (a) total conduc- 
tance, and (b) local charge in the interacting region ^ E^Li ( n i,j)- Here ej/t is taken to be (•) —0.3, (o) —0.7, 
and (A) -1.0. 
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Fig. 10. Contributions of each subband on the conductance of the two-dimensional cluster of the size (M — 4, 
N = 8). Here, e d /t is taken to be (a) -0.3, (b) -0.7, and (c) -1.0. 
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Fig. 11. Contributions of each subband on the local charge in the interacting region J2iLi 2~Hj=i ( n i,j)i °f the two " 
dimensional cluster of the size (M = 4, N = 8). Here, ed/t is taken to be (•) —0.3, (o) —0.7, and (A) —1.0. 




co/t 



Fig. 12. Transmission probability vs ui for the noninteracting system of the size (M = 4, N = 8), for e d /£ = —0.3. 
Here, only the contribution of the mode I (solid line) and mode IV (dashed line) are plotted. The contribution of 
the modes II and III corresponds exactly to the dashed line in Fig. 7 that has plotted for a one-dimensional chain. 
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Fig. 13. Eigenvalues of Tic f° r the two-dimensional Hubbard model of the size (M = 4, N = 8). Here, (a) the 
contribution of the modes I and IV, and (b) that of the degenerate modes II and III. The onsite energy is taken to 
be e d /t = -1.0. 
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Fig. 14. The ed dependence of the conductance through the two-dimensional cluster of the size (M = 4, TV = 8). 
The upper panel (a) shows the total conductance, where the Hartree-Fock result (dashed line) is also plotted for 
comparison. The lower panel (b) shows the contribution of each of the modes. The Coulomb interaction is taken 
to be U/(2nt) — 0.5, so that the electron-hole symmetry holds at ed = ~U/2 (~ — 1.6 i). 



